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Abstract 

We obtained a wide-band spectrum of the Compton-thick Seyfert 2 galaxy Mrk 
3 with Suzaku. The observed spectrum was clearly resolved into weak, soft power- 
law emission, a heavily absorbed power-law component, cold reflection, and many 
emission lines. The heavily absorbed component, absorbed by gas with a column 
density of 1.1x10^^ cm~^, has an intrinsic 2-10 keV luminosity of ~1.6xl0'^^ erg 
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s~^, and is considered to be direct emission from the Mrk 3 nucleus. The reflection 
component was interpreted as reflection of the direct light off cold, thick material; the 
reflection fraction R was 1.36±0.20. The cold material is inferred to be located > 1 
pc from the central black hole of Mrk 3 due to the low ionization parameter of iron 
< 1 erg cm s~^) and the narrow iron line width [a < 22 eV). A Compton shoulder to 
the iron line was detected, but the intensity of the shoulder component was less than 
that expected from spherically distributed Compton-thick material. The weak, soft 
power-law emission is considered to be scattered light by ionized gas. The existence 
of many highly-ionized lines from O, Ne, Mg, Si, S, and Fe in the observed spectrum 
indicates that the ionized gas has a broad ionized structure, with ^=10-1000. The 
scattering fraction with respect to the direct light was estimated to be 0.9±0.2%, 
which indicates that the column density of the scattering region is about 3.6 x 10^^ 
cm~^. This high-quality spectrum obtained by Suzaku can be considered a template 
for studies of Seyfert 2 galaxies. 

Key words: galaxies:active — galaxies: Seyfert — galaxies: individual (Mrk 3) 

1. Introduction 

Active galactic nuclei emit huge amounts of energy resulting from accretion of matter 
onto supermassive black holes. Strong hard X-ray emission from these nuclei is thought to 
be non-thermal, and represented by a power-law with a canonical photon index of ~1.9. The 
strong emission irradiates material around the supermassive black hole, leading to reprocessing 
of the emission. Observing the reprocessed emission can reveal the structure of the nucleus, 
and indicate distribution, abundances, and ionization state of the circumnuclear material. One 
problem is that the reprocessed emission can easily be dominated by the strong nuclear emission, 
e.g., in Seyfert 1 nuclei. The X-ray emission from Seyfert 2 nuclei, however, is known to be 
attenuated by thick matter (e.g. Awaki et al. 1991). Therefore, Seyfert 2 galaxies are among 
the best targets to explore the reprocessing environment around the central black hole. 

Mrk 3 is an X-ray bright Seyfert 2 galaxy at a redshift z=0.013509 (Tifft, Cocke 1988). 
Ginga detected that direct light from the Seyfert 2 nucleus was heavily obscured by thick ma- 
terial with a column density (A^h) of 7x10^^ cm~^ (Awaki et al. 1990). Its intrinsic luminosity 
in the 2-10 keV band was estimated to be 2x10^^ erg s~^, which was comparable to those of 
Seyfert 1 galaxies. Due to the obscuration of the direct light, soft X-ray emission with many 
highly-ionized lines was detected (e.g. Iwasawa et al. 1994, Sako et al. 2000). A Chandra 
observation found that the soft X-ray emission was extended along the [OIII] ionization cone, 
and found that the ionized lines had signatures of emission from a warm, photoionized region 
(Sako et al. 2000). The obscuration of the direct hght has also allowed us to detect a reflec- 
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tion component in the hard X-rays of Mrk 3 (e.g. Turner et al. 1997, Cappi et al. 1999). 
The observed X-ray properties are consistent with Seyfert unification schemes (e.g. Antonucci, 
Miller 1985). From these previous observations, it is known that Mrk 3 has a complex spec- 
trum, which consists of a weak power-law emission, heavily-absorbed power-law emission, cold 
refiection, and many emission lines, but it is difficult to resolve the X-ray spectrum into each 
component. BeppoSAX demonstrated the power of broad-band observations by disentangling 
three of the components in Mrk 3: the soft X-ray emission, the heavily-absorbed component, 
and cold refiection accompanied by an iron line (Cappi et al. 1999). High-quality spectra of 
Mrk 3 were obtained with XMM-Newton, and ionized gas and cold, thick material around the 
nucleus were investigated through measurements of highly-ionized lines and Fe K lines, respec- 
tively. A Compton shoulder to the iron line due to Compton down-scattering gives us evidence 
of refiection by Compton-thick material (e.g. Matt 2002, Watanabe et al. 2003). The intensity 
ratio between Fe-Ka and K-(3 lines is a powerful tool to measure the ionization state of iron in 
the refiecting matter. The observations by XMM-Newton indicated that iron line diagnostics 
have the potential to investigate the cold refiecting matter (Pounds, Page 2005; Bianchi et al. 
2005). 

In this paper, we analyze a wide-band spectrum of Mrk 3 obtained by Suzaku (Mitsuda 
et al. 2007), with the aims of resolving the spectrum into each component and estimating 
properties of the matter around the nucleus. Errors correspond to the 90% confidence level for 
one interesting parameter unless we explicitly note a different confidence level. 

2. Observation and Data Reduction 

Mrk 3 was observed with the Suzaku satellite on 2005 October 22-23 during the SWG 
phase. Suzaku has four X-ray telescopes (XRTs; Serlemitsos et al. 2007), with a spatial 
resolution of 2' (half power diameter). The focal plane detectors of the XRTs are X-ray CCD 
cameras, which are referred to as an X-ray imaging spectrometer (XIS; Koyama et al. 2007). 
Three of the four XISs (XISO, 2, and 3) are front-illuminated CCDs (FI) and the other XIS 
(XISl) is a back-illuminated CCD (BI). The XISs are sensitive to 0.2-12 keV X-rays on a 
18' X 18' field of view. XISl has enhanced soft X-ray sensitivity. Each XIS is equipped with 
two ^^Fe calibration sources, which are illuminated at the CCD corners. Suzaku also has a 
no n- imaging hard X-ray detector (HXD; Takahashi et al. 2007). The HXD has two types 
of detectors, PIN and GSO, and yields sensitivity from ~10 to ~700 keV. Mrk 3 was set on 
the XlS-nominal position, and XIS and HXD data were obtained simultaneously. The XISs 
were operated in the normal clocking mode. The observed XIS and HXD data were processed 
through the pipeline processing, provided as revision 1.2 data. We obtained cleaned event files 
of XIS and HXD, and then analyzed the data following nominal analysis procedures. 
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2.1. XIS Reduction 



The net exposure of the cleaned event data was about 87 ks. First we made X-ray 
images using the cleaned XIS data. Figure 1 shows a 0.4-2 keV band image and a 2-10 keV 
band image. We found seven X-ray sources in the 0.4-2 keV band image. We corrected the 
astrometry of the Suzaku XIS image in adjusting the position of the brightest source in the 
2-10 keV band image on the position of Mrk 3. Table 1 lists the detected sources' positions 
and possible candidates. 

We also list their average count rates, summed over the four XIS detectors in the 0.4-2 
keV band. The observed 0.4-2 keV band X-ray flux of Mrk 3 was estimated to be 5.7x10"^^ 
erg s~^ cm~^. The count rate of the faintest source thus corresponds to 5x10^^^ erg s~^ cm~^ 
assuming that the source has the same spectral shape as Mrk 3 in the 0.4-2 keV band. There 
are bright regions in the corner of the hard band image. These regions were illuminated by the 
^^Fe calibration sources. 

We selected X-ray events within a 4.5' radius centered on Mrk 3 for the spectral analysis. 
There was an X-ray source (IXO 30) in the circular region (e.g. Colbert, Ptak 2002). The 
source was too close to Mrk 3 to separate its X-ray emission from that of Mrk 3. Therefore, 
we included the X-ray events of this source when we extracted the X-ray events of Mrk 3. 
Bianchi et al. (2005) found that the X-ray spectrum of IXO 30 obtained by XMM- Newton was 
represented by a power-law with a photon index of 1.77, and that IXO 30 was non- variable 
over a 2-year monitoring period. Based on their result, the possible contamination from IXO 
30 was estimated to be about 7% of Mrk 3 in the 0.4-2 keV band. In order to estimate an 
X-ray background for Mrk3, we accumulated X-ray events excluding both the source region 
and regions irradiated by the calibration sources (see Figure 1). We obtained both the source 
spectra and background spectra of XIS-0, 1, 2, and 3. 

Response matrices (RMFs) and ancillary response files (ARFs) were generated for each 
XIS independently using xisrmfgen and xissimarf gen in HEASOFT 6.1.2 (Ishisaki et al. 2007). 
The ARF generator takes into account the hydrocarbon contamination on the XIS optical 
blocking filters. To examine the accuracy of the XIS RMFs, we extracted spectra for the 
calibration sources, and then fit the spectra with a two- Gaussian model. The best-fit center 
energies of the calibration sources were about 5 eV lower than the weighted values of the 
theoretical energies for the transitions of Kai and Ka2. The discrepancy of 5 eV was within 
the accuracy (~0.2% at Mn-Ka) of the energy calibration of the XIS. The line widths of the 
calibration sources were found to be (<15) eV. Since the energy scale of each CCD was quite 
similar to within 4 eV, we added the three spectra of the front-illuminated CCDs, XIS-0, 2, and 
3 by using mathpha. Figure 2a shows the summed XIS spectrum and the XIS-1 spectrum, both 
after background subtraction. The XIS team pointed out the existence of an artificial feature 
in the spectrum from 1.825 to 1.843 keV. Thus we removed that energy range from the XIS 
spectra of Mrk 3. In order to fit the summed spectrum, we also added RMFs and ARFs for 
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XIS-0, 2, and 3 spectra into one response file using marfrmf and addrmf in HEASOFT 6.1.2. 
2.2. HXD Reduction 

The HXD-PIN count rate of Mrk 3 was found to be about 0.15 ct s"^, which was about 
30% of the internal (non-X-ray) background (NXB) of the HXD-PIN. PIN background count 
rates are variable and strongly depend on the Suzaku satellite orbit (Kokubun et al. 2007). 
We extracted both the source and background spectra with identical good time intervals. An 
event file of the NXB, reproduced by a model of the NXB for the Mrk 3 observation, was 
provided by the HXD instrument team. The current accuracy of the PIN NXB model for a 1 
day observation is deduced to be about 5% (peak-to-peak residual), which roughly corresponds 
to 2% at the la level. Since the source data was not corrected for instrument dead time, we 
performed the dead time correction by using hxddtcor. The deadtime was 6.5%, and the good 
time exposure came to be about 81 ks. 

Since the HXD-PIN is not an imaging detector, we also had to subtract a contribution 
from the Cosmic X-ray Background (CXB). A response matrix for a fiat field was provided 
from the HXD instrument team. We simulated a CXB spectrum observed by the HXD-PIN 
by assuming the following shape for the CXB: 9.0x10"^ (E/3keV)~°-^^ exp(— i?/40keV) erg 
cm~^ s~^ sr~^ keV~^ (Bolt 1986). Then we added the simulated CXB spectrum to the NXB 
spectrum. Figure 2a shows the PIN spectrum after the background subtraction. The source 
fiux in the 12-70 keV band was 9.4x10"^^ erg cm"^ s~^. 

A response matrix for a point source at the XIS nominal position 
(ae_hxd_pinxinom_20060814.rsp) was used for our spectral analysis. A cross-normalization of 
the HXD-PIN relative to the XIS was derived using Suzaku observations of the Crab. We used 
the cross- normalization of 1.13 at the XIS nominal position (Ishida, the XRT team, 2006). 

We did not use the HXD-GSO data in our analysis, because the background study of 
the HXD-GSO is in progress. When the reproducibility of the GSO background is improved, 
we will include GSO data in order to determine the shape of the X-ray spectrum in the high 
energy band. 

3. Wide-band Spectrum 

A key advantage of the Suzaku observatory is its wide-band spectroscopy from sub-keV 
energies to above 100 keV. We simultaneously fit XIS and HXD spectra with a baseline model 
based on previous observations. The baseline model is defined as follows: 

/(ph s"^ cm-2 keV"^) = e'^^-^^'^'lPLl + e-''^-^^''^PL2 + reflection + ELs], (1) 

where A^hi and A^h2 are column densities for power-law components PLl and PL2, respectively. 
The absorbing cross sections {o'a,i,(Ta,2) are defined in the phabs model in XSPEC vll.3. The 
abundances of Anders & Grevesse (1989) were used. For aa,2, the absorber was placed at the 
cosmological redshift of the source. We fixed the high energy cutoffs of both PLl and PL2 
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at 400 keV, because Cappi et al. (1999) obtained an lower limit to the cut-off energy (i?cut) 
of the power law component E'cut >200keV and our spectral fit could not constrain i?cut- The 
Compton reflection hump was reproduced by pexrav, keeping both the input photon index and 
-E'cut equal to those of PL2. The cosine of inclination and iron abundance in pexrav were fixed 
at 0.45 and one times solar, respectively. The emission lines (ELs) seen in the spectrum were 
represented by the sum of Gaussian components: 

ELs = Y.gauss{E,,auNO, (2) 

i 

where Ei, cri, and Ni are the center energy, width, and intensity of the ith line. For lines below 
2 keV, it is difficult to set a tight constraint on line widths by a Suzaku observation. Thus, we 
fixed them at the values obtained with the Chandra HETG (Sako et al. 2000). The fixed values 
are much smaller than the energy resolution of the XIS (Koyama et al. 2007). The best-fit 
value of A^Hi was less than the Galactic column density 8.7x 10^° cm~^ (Stark et al. 1992), 
thus we fixed the A^hi at the Galactic value. 

The observed spectra were fit with the baseline model in the energy range from 0.4 keV 
to 70 keV, yielding a reduced of 1-15 (d.o.f.=762). Figure 3 shows X-ray spectra in the 
0.4-3 keV and 5.5-7.5 keV bands. Many emission lines are detected. The best-fit parameters 
of the continuum emission components are listed in Table 2, and those of the emission lines are 
summarized in Table 3. An unfolded broad-band spectrum and the baseline model are shown 
in Figure 2b. The baseline model is the same as the best-fit model except only that the line 
width of He-like Fe was the same as that of the Fe Ka line. We included two lines for Si XIII 
and Si XIV in the baseline model, and obtained consistent results with those of XMM-Newton, 
but the result may have some ambiguity due to a uncertainties with the response function 
around the Si edge. We obtained a best-fit photon index and the absorbing column density of 
1.80±0.06 and (1.10±0.06) x 10^^ cm'^, respectively. 

The observed 2-10 keV flux was calculated to be 6.4x10"^^ erg s~^ cm~^, which is 
similar to those of the BeppoSAX (6.5 x 10"^^ erg s"^ cm~^) and XMM-Newton (5.9xl0~^^ erg 
s~^ cm~^) observations, though it is slightly different from the Ginga result ( 8.8x10"^^ erg 
s~^ cm"^). The larger flux in the Ginga observation compared to the BeppoSAX and XMM- 
Newton observations is due to the difference in the absorption column. The intrinsic X-ray 
luminosity of the heavily absorbed component was estimated to be about 1.6x10^^ erg s~^ in 
the 2-10 keV band, assuming Hq = 70 km s~^ Mpc~^ and Aq = 0.73. The X-ray luminosities 
of the other components were estimated to be 1.4x10"^^ erg s~^ in the 2-10 keV band for the 
reflection component and 3.7x10^^^ erg s~^ in the 0.4-2 keV band for the PLl component + 
ELs. 

The reflection ratio R of the cold reflection to PL2 was about 1.36. To estimate the 
error on R, we have to consider the current systematic uncertainties in the background sub- 
traction and the relative PIN/XIS cross-normalization. The statistical uncertainty of the rela- 
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tive PIN/XIS cross-normalisation is only about 3% (Ishida, the XRT team 2006). The change 
in R caused by this uncertainty is smaller than that from the uncertainty of the background 
subtraction. Thus, we evaluated the effect of the NXB reproducibility on R. We modified 
the intensity of the PIN non X-ray background by ±2%, which is about the la level of the 
current reproducibility of the PIN NXB for one day exposure in the 15-40 keV band (Figure 
3 in Mizuno et al. 2006). The net count rates of the source were changed to be +15% and 
-6% for the -2% and +2% NXB, respectively. The best-fit value of R became 1.25 and 1.42, 
respectively. This change is smaller than the statistical error on R, ±0.20, determined by the 
fitting procedure. The value of R derived from the Suzaku observation is different from the 
BeppoSAX value {R = 0.95lJ]i^). 

To estimate the scattering fraction, we have to subtract the contamination from IXO 30. 
Since the 2-10 keV flux of IXO 30, as observed by XMM-Newton, was ~8.5xl0-^'^ erg cm'^ 
s~^ (Bianchi et al. 2005), the contamination of IXO 30 into PLl was estimated to be about 
17%. Thus the ratio of PLl with respect to PL2 was found to be 0.9±0.2%. 

We compared the center energies of ELs with those obtained by XMM-Newton and 
Chandra, and found that the center energies were consistent with those values within 20 eV. 
We also found that our line intensities were more comparable to those for XMM-Newton rather 
than Chandra. Figure 4 shows the comparison between our result and the Chandra result. 
The difference is larger in the soft band and becomes unity above 2 keV. Sako et al. (2000) 
pointed out the spatial extent of the soft band image, which laid approximately in the cross 
dispersion direction of the HETGS. Since Sako et al. (2000) applied a 8 pixel filter in order to 
reduce the noise level in the dispersed spectrum, their spectrum did not cover all the extended 
emission. The consistency with the XMM-Newton results can be explained by the fact that the 
line-emitting region is extended. The equivalent widths of these lines to the PLl and to the 
reflection are listed in Table 4. For the iron line, the line width of about 20 eV (< 29 eV) and 
the intensity of about 4.7 xlO~^ ph s~^ cm~^ are more consistent with the Chandra grating 
results rather than the XMM-Newton result (Bianchi et al. 2005). The equivalent width of the 
iron line to the reflection component was estimated to be ~855 eV. 

3.1. A Compton Shoulder (CS) 

A Compton shoulder component is often represented by a second Gaussian component 
with a width of ~80 eV at an energy about 0.1 keV lower than the primary Fe Ka Gaussian 
(e.g. Matt 2002). However, we represented the component by a pulse function with a pulse 
width of 156 eV in this paper, since the shape of the shoulder component is more similar to 
a pulse function rather than a Gaussian function (e.g. Matt 2002, Watanabe et al. 2003). 
Furthermore, the pulse width of 156 eV corresponds to the maximum energy shift due to a 
single Compton scattering. By fixing the pulse width at 156 eV in our fitting procedure, the 
CS can be represented by only one parameter, /, which is the fraction of the CS intensity with 
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respect to that of the primary Gaussian component. We fit the iron hne bandpass of Mrk 3 
with the Gaussian plus pulse components, finding that the pulse component was detected at 
the 90% confidence level. In this fitting procedure, both the parameters of many lines except 
of the iron Ka line and the photon index were fixed at the best-fit values. The best-fit ratio / 
of the pulse component to the Gaussian component was 10±8%. 

When including the CS, the best-fit values of the line center energy and the line width 
moved slightly to 6.406 keV and (<22) eV, respectively. The line width of cr < 22 eV 
corresponds to a velocity of < 1000 km s~^ due to Doppler broadening. The line width of 20 
eV in the baseline model can be explained by the existence of the CS. The effects of the CS to 
the best-fit center energy and width are shown in a contour map (Figure 5). 

3.2. Iron Ka and KjS Lines 

By determining the center energies and intensities of the iron Ka and K/9 lines accurately, 
we can measure the ionization state of the line-emitting gas (e.g. Bianchi et al. 2005, Yaqoob 
et al. 2007). At first, we obtained a line intensity ratio iRp/lKa of 0.10±0.03 in the baseline 
model (no CS model). The ratio was determined more accurately than obtained by Bianchi et 
al. (2005). In comparison with a calculation of the line ratio in the low ionization region of Fe 
by Palmeri et al. (2003), we constrained the ionization state to Fe I-IX. Figure 6 shows the 
confidence contours of the center energies of Ka versus K^. The calibration uncertainty of the 
XIS response at 6 keV is about 15 eV (Koyama et al. 2007). The center energies of the iron 
Ka and K/5 lines may be shifted within the calibration uncertainty. We shifted the contour 
map by ±7.5 eV. If the contour maps is shifted by -3 eV, the filled circles of Fe II to Fe X 
are included in the 90% confidence region. The line ratio and the line energies indicate a low 
ionization state of the reflection matter. As noted above, the central energy of the iron Ka line 
shifted from 6.401 keV to 6.406 keV when the CS was included into the iron Ka profile. 

3.3. Advanced Analysis of the Reflection Component 

A Compton reflection component with a reflection factor i? of ~ 1.36 was clearly detected. 
Although we assumed one solar abundance of iron in the baseline model, we can measure the iron 
abundance (^peO from the edge seen in the reflection spectrum. We made A^^^ of pexrav free in 
our fitting procedure, and obtained the best-fit value of Ap^ = 0.56±0.10 with an improvement 
of Ax^ = -32. Because an edge structure was also produced by the heavy absorption with 
Nu2, the best-fit Ag was coupled with that (^Fe*^) of the absorption matter. The dependence 
of these abundances is shown in a confidence map (see Figure 7). If the reflecting matter is 
identical to the absorbing matter, A^p^ should be similar to be A^^. Therefore, we estimated 
the iron abundance assuming Af^ = Afj". The abundance was found to be 0.84±0.05 with 
Ax' = -28. 

We measured the energy of the Fe edge seen in the reflection component. Since there 
is no parameter on the edge energy in pexrav, we changed the parameter redshift in the pexrav 
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model in order to obtain the confidence region of tfie edge energy. Tlie best-fit value of the 
redshift was 0.0141;o;qo9, which corresponds to a confidence region of the edge energy of 7.1lto;o5 
keV after correcting for the cosmological redshift. 

We examined the possibility of obscuration of the reflection component in the baseline 
model, because if the covering factor of the reflecting matter is large, the reflection component 
may be obscured by itself. Including the absorption N^^ of the reflection component in the 
baseline model, the fit was significantly improved, with Ax^ ~ -7. The -F-test probability for 
the improvement of is estimated to be about 1%. The best-fit photon index changed to 1.83 
in the case of A^^'^^ =1.7x10^^ cm~^. If the iron abundance was changed to be subsolar instead 
of one times solar, obscuration of the reflection component was not required (A^f^ < 0.9 xlO^^ 
cm~^). 

3.4- A Highly-Ionized Iron Line 

A highly-ionized iron line at the energy of 6.682lto;o23 keV is detected at the ~3cr level. 
This energy is consistent with the energy of a resonance {E = 6.700 keV) or an intercombination 
line {E = 6.675 keV) rather than a forbidden line [E = 6.637 keV). The intensity of this line 
was found to be 3.9l};2xl0~^ ph cm~^ s~^. The observed equivalent width with respect to the 
scattered light (PLl) was about 700 eV. Netzer et al. (1998) estimated the equivalent width 
of a highly-ionized iron line to the scattered light. Our observed equivalent width is consistent 
with their estimation of 500-1000 eV. 

We examined the existence of the hydrogen-like iron line at E = 6.966 keV in the rest 
frame, but we obtained only an upper limit of 2x10"^ ph cm~^ s~^. 

3.5. Absorption model for Nu > ICF^ cm~'^ 

The PL2 source is heavily obscured by thick matter with A^h2=1-1 x 10^^ cm~^. For A^h 
> 10^^ cm~^, the effect of the Compton scattering have to be considered, because optical depth 
of the Compton scattering is close to unity (Ghisellini et al. 1994, Yaqoob 1997). Since the 
phabs in XSPEC considers only photo-electric absorption, we tried to use plcabs in XSPEC 
which describes X-ray transmission of an isotropic source of photons located at the center of 
a uniform, spherical distribution of matter, correctly taking into account Compton scattering. 
This model is not valid above 18 keV, but it is useful to know the effect of the Compton 
scattering. 

The best-fit parameters of the continuum components are listed in Table 5. The fit 
with plcabs is statistically better than the our baseline model (A^^ ~ —10). This model gives 
the best fit column density A^h2 of 1.1x10^^ cm~^ which is consistent with that obtained with 
the baseline model, and a steep photon index of 1.88. The reflection ratio and the scattering 
ratio are changed to i?~l.l and Cf~0.6% due to the increase of the intensity of the direct 
component. The intrinsic luminosity strongly depends on the best fit model. The 2-10 keV 
intrinsic luminosity with plcabs model is estimated to be 2.2x10^'^ erg s~^, which is larger than 
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that obtained by the basehne modeL 

For the emission hnes, most of the best parameters do not significantly change. The 
best-fit center energy of Fe Kf3 line reduced to be 7.04±0.03 with the intensity of 3.9x10"^ ph 
cm~^ s~^. The energy of the Kf3 line is still consistent with that for neutral Fe. 

4. Timing Analysis 

Thanks to Suzaku's wide bandpass, we detected a direct emission component from Mrk 
3. We made a light curve in the 15-40 keV band with the Suzaku PIN to examine the time 
variability of the direct component, which is dominant in this energy band. Background sub- 
traction and deadtime correction were performed for the light curve. The light curve is shown 
in Figure 8. The systematic error of the non X-ray background in the 15-40 keV band for a 
5760 s bin was about 6% (Icr) of the non X-ray background (Mizuno et al. 2006). Considering 
this systematic error, we found that the light curve was stable during the observations, with 
a reduced of 0-40 (d.o.f.=32) against a constant model. The 90% upper limit to the time 
variability was estimated to be cr^j^g<0.01. 

We examined the long-term variability of Mrk 3 by comparing the Suzaku result with 
the BeppoSAX result, obtained during 1997 April 16-18. The 2-10 keV intrinsic luminosity 
observed by BeppoSAX can be estimated from the 0.1-150 keV luminosity by Cappi et al. 
(1999). Although the derived luminosity depends on the best fit model (e.g., it depends on 
values of r, A^h, and R), we obtained an X-ray luminosity of 2.7x10^^ erg s~^. We also compared 
the Ginga and XMM-Newton results, although these results have some ambiguity due to being 
narrow band observations. The Ginga luminosity during the period of 1989 September 27-28 
and the XMM-Newton luminosity during 2000 October 19-20 were estimated to be 2.3x10"^^ 
erg s~^ (e.g. Awaki et al. 1991, Griffiths et al. 1998) and 2.2x10"^^ erg s"^ (e.g. Bianchi et 
al. 2005), respectively. Long-term variability of Mrk 3 was pointed out by Cappi et al. (1999). 
There may be a gradual decrease in the intrinsic luminosity between the BeppoSAX, XMM- 
Newton and Suzaku observations. Further observations with Suzaku are crucial to quantifying 
the long term variability of the intrinsic luminosity of Mrk 3. 

5. Discussion 

5.1. Continuum Emission 

We obtained a wide-band spectrum covering energies from 0.4 to 70 keV. The contin- 
uum emission was resolved into three components: a soft, weak power-law, a cold- reflection 
component, and a heavily absorbed power-law component. We reproduced well the spectrum 
assuming the same photon index for each component. 

The heavy absorbed component with Ah2 ~ 1.1x10^^ cm~^ is considered to be a direct 
component from the nucleus, since the photon index of ~1.8 is similar to the canonical value 
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for Seyfert 1 galaxies and the intrinsic X-ray luminosity in the 2-10 keV band was estimated 
to be 1.6x10^^ erg s~^. The column density we obtained is consistent with the BeppoSAX and 
XMM-Newton observations: For the last decade, Mrk 3 has remained a Compton-thick object. 
Some time variability of the absorbed component is expected, similar to Seyfert 1 galaxies 
(e.g. Markowitz et al. 2003). Mrk 3's intrinsic luminosity has evolved somewhat over the 
last decade. However, short-term time variability was not significantly detected over the 2 day 
observation. The black hole mass (Mbh) has been estimated to be 4.5x10^ Mq from optical 
observations (Woo, Urry 2002). a^^^ for a 2 day observation was expected to be ~3xl0~^ from 
the relation between Mbh and a'^^^ (Papadakis 2004); this is consistent with the upper limit 
obtained during the Suzaku observation. 

5.2. Cold Reflection 

The cold reflection component was detected, with a reflection factor i?=1.36±0.20, and 
the observed X-ray flux in the 2-10 keV band was (3.0±0.1)xl0~^^ erg s~^ cm~^ in the baseline 
model. The value of R is different from that from the BeppoSAX observation, 0.95lo;i5, and 
the observed flux in the 2-10 keV band is consistent with the X-ray flux of ~ 3.4x10^^^ erg 
cm~^ s~^ measured during the XMM-Newton observation (Bianchi et al. 2005). The direct 
nuclear emission may be gradually decreasing from the time of the BeppoSAX observation, to 
the XMM-Newton and Suzaku observations. If the reflector is a pc scale torus, then a time 
delay between the reflection component and the direct component is likely to be on a timescale 
of years. Continued X-ray spectral monitoring of Mrk 3 with Suzaku would thus be useful for 
constraining the location of the reprocessing matter. 

We detected the CS of the iron Ka line with a 90% confidence level. It is known that 
the intensity ratio / depends on the column density of the scattering/reflecting material (e.g. 
Matt 2002, Watanabe et al. 2003). The value of / = 0.1 roughly corresponds to an optical 
depth of Compton scattering: t ^ N^, where ctt and are the Thomson cross-section 
and the column density of the reflecting matter, respectively. was estimated to be 2 xlO^^ 
cm~^, which is smaller than the absorbing column A^h2- This means that a uniform density 
spherical distribution of reflecting matter around an X-ray emitter was not appropriate. Matt 
(2002) calculated the CS properties for an isotropically illuminated plane-parallel slab. In this 
geometry, / is smaller than that in a spherical distribution. In order to explain a small value 
of /, we have to consider a geometric configuration wherein material has a large optical depth 
before scattering or absorption by iron, and a small optical depth for reprocessed light. A torus 
geometry is probably worth studying for at least part of the Compton shoulder emission. 

The ionization parameter, ^, of the reflecting matter can be estimated from the iron Ka 
and Kj3 lines. From the observed intensity ratio of these lines, a low ionization state, Fe I-IX, 
has already been pointed out. We examined the low ionization state using the center energies of 
these lines. Palmeri et al. (2003) calculated theoretical energies of the Fe Ka and K/5 lines for 
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Feii-Feix- We plotted the theoretical energies on Figure 6. The plot also demonstrates that the 
ionization state of iron is lower than Fe IX. The low ionization state of iron was also confirmed 
with the iron edge energy. Kallman, Baustina (2001) and Kallman et al. (2004) estimated ion 
fraction of a low-density ionized gas as a function of ionization parameter C, ■ Their estimate 
is valid for A^h < 10^^ cm~^. We suppose that the absorber is the same as the reflector. In the 
case of a pc scale torus, the mean density of the torus is estimated to be 3x10^ cm~^. The 
ionization parameter with the largest emissivity for Fe IX is estimated to be about 1 erg cm 
s~^. The luminosity (L'°°) to be used when estimating the ionization parameter is defined as 
the luminosity in the 1 - 1000 ryd band. Assuming a simple power-law model with an energy 
index of 0.8, L^""^ was deduced to be 3 times the luminosity in the 2-10 keV band. From the 
definition of ^, the distance r from an X-ray source with L'°° was deduced as follows: 

r ion r ion „ p 

r = _ 2 V ^1/2/ ^ ^~l/2( ^ )-^/^(nc) r3l 

■ ^5x 1043 erg s-1^ 4o6 cm-3^ ^ergcms^i^ ^ 

The ionization parameter indicates that the reflecting region is located at r > 1 pc from the 

X-ray emitting region. 

The location of the reflecting region was also estimated from the iron line width, assuming 

Keplerian motion of the iron line-emitting matter and Doppler broadening of the iron line. The 

velocity (vfwhm) estimated from the line width is expected to be about 2fksinz, where is 

the Keplerian velocity and i is the angle between the rotation axis and our line of sight. The 

distance r from the center of gravity is deduced as follows: 

^^^ 1300 kms-^ ^ M3smz(pc), (4) 
where Mg is the total amount of mass within r in units of 10® M©. Assuming that the central 
black hole mass dominates the total mass, the central black hole mass of 4.5x10® Mq and the 
upper limit to the line width of 2400 km s~^ yield an estimate on the distance: r > 1.3 sin^z 
pc. Assuming that the inclination angle i is equal to the narrow line region (NLR) inclination 
angle of 80-85° (Ruiz et al. 2001), the distance is greater than 1.3 pc, which corresponds to 
~3xl04 Ts, where Ts is the Schwarzschild radius for a black hole mass of 4.5x10® Mq. The 
estimate using the iron line width also suggests that the reflecting region is far from the central 
black hole. 

The iron abundance of the reflecting region was estimated from the iron edge structure 
seen in the reflected component, and an abundance of ~0.56 was obtained for cost = 0.45. 
Linking it with the abundance of the absorbing column density A^h2, the abundance changed 
to ~0.84 solar. Metal abundances were also estimated from the equivalent width of the iron 
Ka line with respect to the reflection component. For cosz=0.45 and a metal abundance of one 
solar, the equivalent width of an Fe Ka line with respect to a reflection component is derived 
to be about 1200 eV assuming a plane-parallel, semi-inflnite slab of neutral matter (e.g.. Matt 
et al. 1997, Pounds & Page 2005). Matt et al. (1997) used the Anders & Grevesse abundances 
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(hereafter AG abundances), which is the same as what we used in XSPEC. The observed 
equivalent width of 855±35 eV corresponds to an iron abundance of 0.7±0.1 solar, which is 
consistent with the iron abundance estimated from the edge structure. The iron abundance in 
Mrk 3 is thus estimated to be sub-solar, using the AG abundances and the plane-parallel slab 
geometry. We compared the abundance with those of other AGNs observed by Suzaku. The 
iron abundances in MCG-5-23-16 and Gen A were measured to be 0.65±0.15 solar (Reeves 
et al. 2007) and 1.2±0.1 solar (Markowitz et al. 2007), respectively. However, since those 
authors used different abundance tables (Wilms, Allen, & McGray (2000) and Lodders (2003), 
respectively), we converted them to the iron abundances in the case of the AG abundances. 
The converted abundances were ~0.4 and ~0.8 solar in MGG-5-23-16 and Gen A, respectively. 
The iron abundance in Mrk 3 is therefore similar to those of these AGNs; taking the differences 
between various abundance table into account, the iron abundance of Mrk 3 ranges from 0.5 to 
~1 solar. 

We estimated metal abundances of other elements from equivalent widths of those lines 
with respect to the reflection component. Matt et al. (1997) also calculated the equivalent 
widths of Si and S Ka lines to be 240 eV and 165 eV for cosz=0.45 and a metal abundance 
of one solar, respectively. Gomparing these predictions to the observed equivalent widths of 
300±35 eV and 160±30 eV, respectively, for these lines, the metal abundances of Si and S are 
estimated to be about one solar. We note that there is some ambiguity for the observed intensity 
of the Si Ka line, but the intensity is consistent with that obtained using XMM-Newton (see 
Table 3). If the Si and S Ka lines are emitted at the reflection region along with the Fe lines, the 
abundance ratios of Z^i/Zpc and Z^/Zp^ are 1.7±0.2 and 1.3±0.2, respectively. These results 
are similar to those obtained for Gen A by Markowitz et al. (2007). The abundance ratios 
suggest that the metals may be produced by type II SNe (e.g. Tsuru et al. 2007). However, 
we have to consider the possibility of contamination of the Si K and S K lines from emission 
due to a photoionized region with a low ionization parameter, since the photoionized region 
has a broad ionization structure, as discussed below. This possibility could be related to the 
difference in observed Si Ka intensities between the Suzaku and Ghandra observations. 

5.3. The Scattering Region 

The soft continuum emission, described by the power-law component PLl, is considered 
to be scattered light from the photoionized region, as with polarized light in the optical band. 
Thanks to Suzaku's wide bandpass, the ratio C/ of the scattered to direct components is 
accurately estimated to be 0.9±0.2%. This ratio is consistent with that assumed by Miller & 
Goodrich (1990). The ratio approximately corresponds to r AQ/Arc, where AQ is the covering 
factor of the scattering region, and r is the scattering optical depth. The column density (A^^*^*) 
for the scattering is described as r/ax, where ax is Thomson cross-section. 

^r = ^ = 1.3xlO-(^)(^)-W^ (5) 
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The opening angle of the scattering region may be estimated from the opening angle of the NLR. 
Capetti et al. (1995) found a NLR opening angle of >100°. Assuming a half opening angle 50° 
for the bipolar ionization cone (Capetti et al. 1995), N^^^ is deduced to be 3.6x10^^ cm~^. The 
half opening angle of 50° is larger than the estimate by Ruiz et al. (2001) due to the inclusion 
of all the Z-shape emission components in the NLR in the estimates of Capetti et al. (1995). 
PLl had an identical photon index to PL2, and had a low absorption column corresponding 
to A^H ~ 6.4x10^° cm~^. This suggests that the scattering region is highly ionized. In the 
case of a uniform electron density, the size of the scattering region is estimated to be about 
50 (L/5xl0^^ erg s~^) (Ah/3. 6x10^^ cm"^)"^ (,^/10)"^ pc, and the density in the scattering 
region is estimated to be about 240 cm~'^. 

Information on the ionized region can be obtained from the emission-line spectrum. A 
photoionized plasma can be characterized by an ionization parameter ^. Since He-like iron and 
He-like oxygen do not coexist in a photoionized plasma with a uniform ^ (Kallman, Bautista 
2001), the ionized gas in Mrk 3 must have a broad ionization structure with (^=10-1000, as 
pointed out by previous observations (e.g. Pounds, Page 2005). 

As mentioned above, the scattering region may extend up to 1 kpc (Sako et al. 2000). 
This is inconsistent with the estimate of the uniform electron density model. We tried to 
explain the extent of the scattering region by assuming that the density gradually decreases as 
the distance from the nucleus increases, i.e. n=nQr~^. The column density was derived from 
the integral of the density from r=r[^ to r^nt, where and r^nt are the inner and outer radii of 
the scattering region, respectively. In r < 1 pc, the scattering region might be obscured by the 
dusty torus. We thus set rin=lpc. For rout, "we adopted the extent observed by Chandra. For 
Ah=3.6x10^^ cm~^, the electron densities at 1 pc and 1 kpc were estimated to be 1700 cm~^ 
and 1.7 cm~^, respectively. The ionization parameters ,^ at 1 pc and 1 kpc are 3300 and 3.3, 
respectively. This is consistent with the broad ionization structure. 

6. Summary and Conclusion 

We obtained a wide-band spectrum of Mrk 3 with Suzaku. The spectrum in the 0.4-70 
keV band was described by our baseline model, which consists of a power-law emission com- 
ponent, heavily obscured power-law emission, a cold reflection component, and many emission 
lines. The heavily absorbed component is considered to be a direct nuclear emission compo- 
nent. The photon index of the obscured component was 1.80±0.06. The intrinsic luminosity 
of the direct component was estimated to be 1.6x10^^ erg s~^ in the 2-10 keV band. We did 
not find any significant time variability during the 2-day observation, but there is evidence 
for long-term variability from a comparison with results of previous observations. The nuclear 
activity in the 2-10 keV band is almost completely obscured by thick matter with Ah2 ~1.1x 
10^^ cm~^, which suggests that Mrk 3 was in a Compton-thick phase. 

The cold reflection component was detected with a reflection factor of R ~1.36. The 
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observed luminosity was estimated to be 1.4x10^^ erg s~^ in the 2-10 keV band. A strong iron 
Ka line was detected with the EW of ~855 eV with respect to the cold reflection. We examined 
the existence of a Compton shoulder to the iron Ka line, and detected the shoulder with 90% 
confidence. However, the intensity of the shoulder component was smaller than that expected 
from Compton reflection by spherically distributed thick matter. To resolve this issue, we need 
to consider a suitable distribution for the thick matter. We evaluated an ionization parameter 
^ of the reflection component by analyzing iron Ka and K(3 lines in detail. The ionization state 
of iron was determined to be lower than Fe IX. This was also confirmed from the energy of the 
iron edge in the refiected component. The low ionization state can be realized for ^ < 1 in a 
photoionized plasma. The small ^ suggests that iron emitting matter, considered to be refiected 
matter, is located at > 1 pc from the strong X-ray emitter. This result is also deduced from 
the iron line width. The large distance are consistent with a dusty torus. The iron abundance 
was estimated by measuring the depth of the iron edge structure in the refiected component 
and by evaluating equivalent widths of a Ka line from cold matter. Both results show that 
iron abundance is sub-solar assuming the Anders & Grevesse abundances. Suzaku satellite 
can determine the EWs of emission lines with respect to the refiection component due to the 
Suzaku wide bandpass accurately. 

We detected a weak power-law emission component PLl and many highly-ionized lines, 
which are associated with a highly-ionized plasma which is photoionized by the inferred strong 
nuclear continuum. The weak power-law emission component was considered to be scattered 
light, with a scattering efficiency estimated to be about 1%. The efficiency corresponds to a 
value for the column density of the ionized matter of about 3.6 xlO^^ cm~^. The ionized gas has 
a broad ionization structure, ^ = 10-1000, deduced from the detection of both He-like oxygen 
and He-like iron lines. The size of the scattering region was deduced to be about 50 pc, using a 
simple uniform density model. However, to explain the spatial extent of the scattering region, 
the density may gradually decrease with as the distance from the nucleus increases. 

The authors wish to thank the staff of the Suzaku team. We also thank an anonymous 
referee for useful comments. This study is carried out in part by the Grant-in-Aid for Scientific 
Research (17030007 H.A. and 17740124 Y.T.) of the Ministry of Education, Culture, Sports, 
Science and Technology. 



15 



References 



Anders, E. &; Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Antonucci, R. & Miller, J.S. 1985, ApJ, 297, 621 

Awaki, H., Koyama, K., Kunieda, K., & Tawara, Y. 1990, Nature, 346, 544 
Awaki, H., Koyama, K., Inoue, H., & Halpern, J.P. 1991, PASJ, 43, 195 
Bianchi, S., Minuitti, G. Fabian, A.C., & Iwasawa, K. 2005, MNRAS, 360, 380 
Bolt, E. 1987, lAU Symp., 124, 611 

Capetti, A., Macchetto, F., Axon, D.J., Sparks, W.B. & Boksenberg, A. 1995, ApJ, 448, 600 

Cappi, M. et al. 1999, A&A, 344, 857 

Corbert, E.J.M.. & Ptak, A.F. 2002, ApJS, 143, 25 

Ghisellini, G., Haardt, F., Matt, G., 1994, MNRAS, 267, 743 

Griffiths, R.G., Warwick, R.S., Georgantopolos, I., Done, C., & Smith, D.A. 1998, MNRAS, 298, 1159 
Ishida, M., & the XRT team, 2006, JX-ISAS-SUZAKU-MEMO-2006-40, "Relative normalizations of 

XIS and HXD-PIN for Vl.2.2.3 data" 
Ishisaki, Y. et al. 2007, PASJ, 59S, 117 

Iwasawa, K., Yaqoob, T., Awaki, H., & Ogasaka, Y. 1994, PASJ, 46, L167 

Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 

Kokubun, M. et al. 2007, PASJ, 59, S53 

Koyama, K. et al. 2007, PASJ, 59, S23 

Lodders, K. 2003, ApJ, 591, 1220 

Markowitz, A. et al. 2003, ApJ, 593, 96 

Markowitz, A. et al. 2007, astro-ph/07043743 

Matt, G., Fabian, A.C. & Reynolds, C.S. 1997, MNRAS, 289, 175 

Matt, G. 2002, MNRAS, 337, 147 

Miller, J.S., & Goodrich, R.W. 1990, ApJ, 355, 456 

Mitsuda, K. et al. 2007, PASJ, 59, SI 

Mizuno, T. et al. 2006, JX-ISAS-SUZAKU-MEMO-2006-42, "Reproducibility of the HXD PIN Non 

X-ray background" 
Netzer, H., Turner, T.J., & George, I.M. 1998, ApJ, 504, 680 
Palmeri, P., Mendoza, C., Kallman, T.R., & Bautista, M.A. 2003, A&A, 410, 359 
Papadakis, I.E. 2004, MNRAS, 348, 207 
Pounds, K.A., & Page, K.L. 2005, MNRAS, 360, 1123 
Reeves, J.N. et al. 2007, PASJ, 59, S301 

Ruiz, J.R., Crenshaw, D.M., Kraemer, S.B., Bower, G.A., Gull, T.R., Hutchings, J.B., Kaiser, M.E. 

& Weistrop, D. 2001, AJ, 122, 2961 
Sako, M., Kahn, S.M., Paerels, F., & Liedahl, D.A. 2000, ApJ, 543, L115 
Serlemitsos, P.J. et al. 2007, PASJ, 59, S9 

Stark, A.A., Gammie, C.F., Wilson, R.W., Bally, J., Linke, R., Heiles, C., & Hurwitz, M. 1992, ApJS, 
79, 77 

Takahashi, T. et al. 2007, PASJ, 59, S35 
Tifft, W.G., & Cocke, W.J. 1988, ApJS, 67, 1 

16 



Tsuru, T.G. et al. 2007, PASJ, 59, S269 

Turner, T.J., George, I.M., Nandra, K., & Mushotzky, R.F. 1997, ApJS, 113, 23 

Watanabe, S., Sako, M., Ishida, M., Ishisaki, Y., Kahn, S.M., Kohmura, T., Morita, U., Nagase, F., 

Paerels, F., &Takahashi, T. 2003, ApJ, 597, L37 
Wilms, J., Allen, A, & McCray, R. 2000, ApJ, 542, 914 
Woo, J., & Urry, CM. 2002, ApJ, 579, 530 
Yaqoob, T. 1997, ApJ, 479, 184 
Yaqoob, T. et al. 2007, PASJ, 59, S283 



17 



Table 1. X-ray sources detected by Suzaku 



ID Position(J2000) count rate (x 10"^ ct/s/XIS)* possible candidate 



1 


(06'^ 


13™ 


43** 


+71 07 30 ) 


450±1G 


87CjB 060751.7+710809 (I3L LAC) 


2 


(06^ 


13"^ 


46" 


+71° 03' 47") 


59±5 


IWGA J0613.7+7103 (ROSAT X-ray source) 


3 


(06^ 


14m 


53" 


+71° 05' 10") 


5G±5 


2MASX J06145765+7105520 (?) 


4 


(06^ 


15" 


10" 


+71° 07' 34") 


50±4 


UGC 03422 (?) 


5 


(06'^ 


15"" 


35" 


+71° 02' 19") 


2200+24 


Mrk 3 


6 


(06'^ 


16"" 


07" 


+71° 06' 38") 


88+5 


detected by XMM-Newton 


7 


(06'^ 


16"" 


16" 


+70° 56' 06") 


20+4 


IWGA J0616.2+7056 (ROSAT X-ray source) 



*: Average count rates of 4 XIS detectors in the 0.4-2 keV band within a 1''.4 diameter circle after background 
subtraction. 



Table 2. Best-fit model parameters to the continuum emission 



photon index 


APLI 


^ref 


A^H2 


^PL2 id (d.o.f.) 








(cm-2) 




1.80+0.06 


(1.27+0.05) xlO-'^ 


(1.56±0.11)xl0-2 


(1.10+0.06)xl024 


(1.15+0.20)xl0-2 874.8 (762) 



Note — - A^Hi was fixed at the Galactic column density. 
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Table 3. Line energies 



center Giiorj 
(kcV) 



-0.006 
„ ,,,+0.003 
"•''^ -0.004 
„ + 0.005 



c width 
(oV) 

:d) 
:d) 
:d) 
:d) 
:d) 
:d) 
:d) 
:d) 

2 (fixed) 

2 (fixed) 

2 (fixed) 

2 (fixed) 

5 (fixed) 

5 (fixed) 



intensity 
(XIO-"' ph cm-^ 



XMM-Newton RGS.pn^ 
(kaV) 



«-™="or3 

0.562it0.001 

0.57±0.01 
0.655±0.005 
0.741 + 0005 

-o.oos 

O.S27±0.002 
O.i 









N VII Lya (0.5003 keV) 


— 22 


0.561 


37 


OVII Kq f (0.561 koV) 


25+1* 
— 13 

32+6 
— 8 


0.569 
0.654 


<13, <11 
11 


OVII Ka i + r (0.569/0.574 kcV) 
OVIII Lya (0.654 keV) 


28+^, 

— 11 


0.780 


6.3 


O VII RRC ( >0.739 keV) 


11 + ^ 
— 3 


0.824 


4.2 


Fe XVII L (0.826 kcV), NeK? 


16±4 


0.873 


3.9 


OVIII RRC (>0.871) 


18+* 
— 8 


0.905, 0.915 


<2.4, <2.5 


Ng IX Ka f+i (0. 905/0. 915kcV) 


18+^ 
-10 

19+^ 


0.922 
1.022 


6.7 
7.8 


Ne IX Kn r (0.922kcV) 
Ne X Lya(1.022 keV) 


10+8 


1.053 


3.5 


FeXXII L (1.053 kcV) 


8 + 8 
"-4 


1. 170 


1.4 


FgXXIII L (1.170 keV) 


8+^ 


1.211 


2 


Ne X Ly/3 (1.211keV), Mg-K« (1.2536 keV) 


8+^ 


1.332 


<0.8 


Mg XI Ka f(1.331kGV) 


7+7 

-3 


1.352 


2.1 


Mg XI Ka i+r (1.344/1.352 keV) 


R+8 






IVIg XII Lya (1.472 keV) 


7±1 


1.7407 + 0.0073 


2.2 + 0.4 


Si-Ka (1.740 kcV) 


6±1 


1.841 


2.8 


Si XIII f (1.854 kcV) 




1.864 


3.5 


Si XIII r+i (1.865keV) 


2±1 


2.008 


3.5 


Si XIV LyQ(2.005 kcV) 


4+1 






S-Ka (2.307 kcV) 




2.454 


3.7 


S XV(2. 447/2. 461 kcV) 




2.625 


1.8 


S XVI Lya (2.620 kcV) 


38±2 


6.391+0.004 


49 + 5 


Fc-Ka (6.4038 keV) 


4±2 


6.68 


21 


Fe XXV (6.637, 6.675, 6.700 kcV) 


5±2 


7.01 


4.2 


H-likc Fe (6.952/6.973 kcV), Fc K/S 7.058 



0.495^ 



0.566^ 



0.007 



0.659^ 



0.816^ 



0.892^ 



O+0.003 
-0.006 
0.736±0.003 
„+0.004 
-0.005 
0.870±0.009 
, + 0.009 
-0.010 
0.920±0.003 

1.006±0.003 
1.066 + 0.004 
1.152 + 0.006 
, „,„ + 0.008 
''"'-0.006 
1.337+0.004 



1.462 + 0.008 



-0.005 
„ + 0.013 
-0.017 



1.998 + 0.013 

2.315 + 0.011 

2.444+''-02* 
^•***-0.030 
„ „, , + 0.035 
^•''l''-0.040 
6.402+0.003 
, + 0.021 



6.682 
7.090' 
7.424 



0.023 
+ 0.020 
0.022 
+ 0.060 
0.050 



1 (fi, 

1 (fi> 

1 (fi, 

1 (fi, 

2 (fi> 
2 (fi, 
2 (fi, 
2 (fi. 



5 (fi> 
5 (fi> 
5 (fij 



5 (fij 
5 (fi> 



:ed) 
:ed) 
5 (fixed) 
5 (fixed) 
20 (< 29) 
(<42) 
20 



33 + 7 
84 + 20 
52+14 
52 + 5 
54 + 4 
36 + 3 
28 + 7 
31 + 3 
38 + 

26.3 + 2.0 

14.9 + 1.6 

9.8 + 

10.2^ 

11.5"^ 



8 + 5 
-6 



« + 1.4 
-1.2 
, + 1.3 



, + 1.1 
-1.0 



4.5 + 0.9 
7.3 + 0.9 
2.7+ 



,+ 1.1 
-1.0 



2.6 + 0.9 
4.4 + 0.8 
1.4 
1.0 
47.4' 
3.9 
4.9 
3.8 + 1.1 



^+0.8 
-0.9 
„+0.7 
-0.8 
+ 1.8 
2.0 
+ 1.5 
1.2 
+ 1.2 
1.4 



s+0.01 
-0.02 



0.91"' 



,+0.01 
-0.03 
0.924+0.002 
,,„, + 0.05 

i.08+0;0^ 

1.170 



1.24^ 
1.34+ 
1.36 + 
1.47+ 



.j + 0.02 
-0.04 
^+0.01 
-0.04 
R + 0.03 
-0.04 
7+0.01 
-0.02 
1 744+''°l*'t 
-0.002 



1.99;° - t 



6.415+0.006 + 
6.7ll°-03t 

7.06* 
7.60 + 0.05 + 



Ni-Ka (7.4782 kcV) 



*: Corrected for cosmological redshift (z=0.0135) 

f: XMM-Newton EPIC pn 

1 Pounds, Page 2005, 2 Sako et al. 2000 
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Table 4. Equivalent width 



center energy* 


Ewt to PLl 


Ewt to reflection 


ID 


(kcV) 


(cV) 


(i-V) 




-0.006 


69 


670 


N VII Lya (0.5003 kcV) 


551+0.003 
-0.004 


290 


6.7 k 


OVII Ka f (0.561 kcV) 


0.566+°-™S 


120 


2.5 k 


OVII Ka i+r (0.569/0.574 kcV) 


+o:oo3 

-0.006 


190 


2.9 k 


OVIII Lya (0.654 kcV) 


0.736±0.003 


240 


3.4 k 


O VII RRC ( >0.739 kcV) 


„ „,„+0.004 
°'*^''-0.005 


190 


2.0 k 


Fc XVII L (0.826 kcV). NcK? 


0. 870it0.009 




3.0 k 


OVIII RRC (>0 871) 


n „„„+0.009 

"■"^ ''^-0.010 






Nc IX Ka f+i (0 905/0 915kcV) 


0.920±0.003 


350 


3.4 k 


Nc IX Ka r (0.922kGV) 


1.006±0.003 


210 


1.6 k 


Nc X Lya(1.022 kcV) 


1.066±0.004 


130 


900 


FcXXII L (1.053 kcV) 


1.152±0.006 


99 


540 


FcXXIII L (1.170 kcV) 


, „o,+0.008 
^■^■^^-0.006 


130 


580 


Nc X Ly0 (1.211kcV), Mg-Ka (1.2536 keV) 


1.337±0.004 


150 


590 


Mg IX Ka f(1.331kGV) 
Mg IX Ka i + r (1.344/1.352 kcV) 


1.462±0.00S 


68 


210 


]VIg XII Lya (1.472 kcV) 


1 741+0.006 
-0.005 


150 


300 


Si-Ka (1.740 kcV) 


,.,,„+0.013 


62 


120 


Si XIII f (1.854 kcV),Si XIII r+i (1.865kcV) 


1.998±0.013 


68 


100 


Si XIV Lya(2.005 kcV) 


2.315±0.011 


160 


160 


S-Ka (2.307 koV) 


2.444+0024 
^■***-0.030 


53 


47 


S XV(2. 447/2. 461 kcV) 


2.612+0035 
-0.040 


41 


32 


S XVI Lya (2.620 kcV) 


6.402±0.003 


10.8 k 


855 


Fc-Ka (6.4038 kcV) 


+0.021 
-0.023 


960 


68 


Fc XXV (6.637, 6.675, 6.700 kcV) 


7 090+0 020 
'■"""-0.022 


1.30 k 


81 


H-Iikc Fc (6.952/6.973 kcV), Fc K/3 7.058 


7.424+0-060 
-0.050 


1.1 k 


120 


Ni-Ka (7.4782 koV) 



*: corrected for cosmological redshift (z=0.0135) 
f : EW; equivalent width 



Table 5. Best-fit model parameters to the continuum emission with plcabs 



photon index 


ApLl 


•^ref 


Nh2 


ApL2 (d.oi.) 








(cm^^) 




1.88±0.07 


(1.23±0.06)xl0-4 


(1.86±0.10)xl0-2 


(1.12±0.06)xl024 


(1.69±0.17)xl0'2 864.4 (762) 



Note — - A^Hi was fixed at the Galactic column density. 
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#1 



#7 




#3 



1X0 30 



#5 
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Fig. 1. X-ray images of the Mrk 3 region observed by Suzaku. The left and right panels show the images 
in the 0.4-2 keV and 2-10 keV bands, respectively. 




observed energy (keV) 

Fig. 2. Wide-band spectrum of Mrk 3 observed with Suzaku (left). The black, red, and green crosses 
are spectra obtained by the XIS-BI, XIS-FI, and HXD-PIN, respectively. The energy range around the 
Si-K edge (1.875-1.840 keV) is ignored for the spectral fitting. The right figure shows an unfolded X-ray 
spectrum. The yellow, green, and blue lines show the continuum emission: PLl, reflection, and heavily 
absorbed PL2, respectively. The red line shows the best-fit baseline model. 
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observed energy (keV) 

observed energy (keV) 

Fig. 3. X-ray spectra in tlie 0.4-3keV and 5.5-7.5 keV bands. Tlie best-fit center energies of the emission 
lines arc denoted by the green Unes. 
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Fig. 4. Line intensity ratios between Suzaku and previous observations. 
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Fig. 5. The confidence contours of the fraction of Compton shoulder emission with respect to the Ka 
core versus the center energy (left) and the line width (right) of the iron line. The 68%, 90%, and 99% 
confidence levels are plotted. 




E Ka (keV) 

Fig. 6. The 68%, 90%, and 99% confidence contours of the center energies of iron Ka versus K/3. The 
filled circles denote the center energies of these lines for Fe II to Fe X (Palmeri et al. 2003). Since the 
calibration uncertainty of the XIS response at 6 keV is about 15 eV, we show the contour maps with 
energy shifts of +7.5 cV and -7.5 eV. 
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Fig. 7. The confidence contour of iron abundance determined by the edge feature seen in the reflection 
component and the absorption component. The 68%, 90%, and 99% confidence levels are plotted. 
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Fig. 8. The light curve of Mrk 3 in the 15-40 keV band. The error bars show only the statistical error of 
the count rates. A systematic error due to the non X-ray background subtraction is about 0.016 c s^^. 
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